Time dependent spin-dressing using a 3 He atomic beam 
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We have performed high precision experimental measurements of spin precession using a dressed 3 He atomic 
beam. Spin dressing uses an oscillating magnetic field that is both detuned to a high frequency and orthogonal 
to a static magnetic field to effectively change the gyromagnetic ratio of a spin. We verify the validity of 
the spin-dressing Hamiltonian in regions beyond the limiting solution in which the Larmor frequency is much 
smaller than the frequency of the dressing field. We also evaluate the effect of magnetic field misalignment, 
i.e. if the oscillating magnetic field is not orthogonal to the static magnetic field. Modulation of the dressing 
field parameters is also discussed, with a focus on whether such a modulation can be approximated merely as 
a time dependent, dressed gyromagnetic ratio. Furthermore, we discuss implications for a proposed search for 
the neutron electric dipole moment, which would employ spin-dressing to make the effective 3 He and neutron 
magnetic moments equal. 
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I. INTRODUCTION 

The existence of a permanent electric dipole moment 
(EDM) of a particle would be an indication of both parity (P) 
and time reversal (T) asymmetries in the fundamental interac- 
tions that describe the particle. Because of the CPT theorem, 
which states that any Lorentz-invariant theory must conserve 
the combined symmetry operations of charge conjugation (C), 
parity, and time reversal, the existence of an EDM would be 
indicative of CP asymmetry as well [1]. Although cosmolog- 
ical evidence suggests that the fundamental laws of physics 
violate CP symmetry, the standard model of particle physics 
currently does not produce enough CP violation to explain the 
observed matter anti-matter asymmetry in the universe [2]. 
Extensions to the standard model, such as supersymmetry, 
tend to contain more CP violation and therefore produce larger 
EDMs [1]. Such large EDMs could be detected in the next 
generation of experiments. Moreover, experimental searches 
for the neutron EDM (nEDM) place the most stringent lim- 
its on the so-called 6 term in QCD. The lack of this term in 
the standard model Lagrangian has led to the proposal of the 
axion [3-5]. 

A proposed experiment to detect the nEDM would use spin 
polarized 3 He as a comagnetometer [6]. When a neutron is 
absorbed by a 3 He atom, the ensuing nuclear reaction releases 
764 keV of energy, distributed as kinetic energy among the 
charged reaction products (a proton and a triton). If the reac- 
tion occurs in liquid helium, some of the kinetic energy lost 
to the liquid produces scintillation light [7, 8]. Because the 
absorption cross section depends on the relative orientation of 
the two spins [9], i.e. S„ - S3, the scintillation signal in a bath 
of polarized 3 He and neutrons precessing in a magnetic field 
would be time-modulated as cos[(y„ -y^Bt+tp], where y„ and 
73 are the gyromagnetic ratios of the neutron and 3 He, respec- 
tively, and S„ and S3 are the spin vectors of the neutron and 
3 He, respectively. The EDM of the 3 He atom is expected, on 
firm theoretical grounds, to be vastly smaller than the nEDM. 
The presence of an EDM would modify the precession of the 
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FIG. 1. (Color online) Eigenvalues of the spin-dressing Hamiltonian 
(Eq. 2) as a function of dressing parameter x in units of the energy 
of the dressing field photons for y = 0.5. The difference between the 
two central states is associated with the dressed gyromagnetic ratio 
Y for that value of x and y (see text for full explanation). (Inset) 
Appearance of an avoided crossing near x = 2.32 when the angle 
between the dressing field Bj and static field So is 90° (solid line), 
89° (dashed line) and 88° (short dashed line). 



neutrons if an electric field is applied and therefore change the 
modulation frequency of the scintillation signal. 

The magnetic moments of the neutron and 3 He atom are 
equal to within 12% (y„/27r = -2.916 kHz/G [10] and 
ji/2n = -3.243 kHz/G [11, 12], respectively), and so moni- 
toring the differential precession frequency as described above 
reduces by nearly an order of magnitude the undesired effects 
of external magnetic field noise and systematic variations of 
magnetic fields [6]. An ideal experiment would compare a 
neutron to an atom that has exactly the same magnetic mo- 
ment. No such atom exists, and furthermore the only atom that 
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FIG. 2. (Color online) Effective gyromagnetic ratio as a function 
of the dressing field amplitude, x. Each curve represents a different 
value of y; the values of y that are plotted are 0, 0.2, 0.4, 0.6, and 
0.8. The arrow indicates the direction of increasing y. The curve 
for y = is merely Jq(x), which is the analytic solution to the spin- 
dressing Hamiltonian. 



will remain in solution with liquid helium at low temperatures 
is 3 He. However, it is possible to tune the effective magnetic 
moment of a particle by the use of spin-dressing, and, more- 
over, to make the effective gyromagnetic ratios of the neutron 
and 3 He atom equal. This procedure is referred to as critical 
dressing. The parameters to achieve critical dressing depend 
only weakly (to the second order) on the static magnetic field. 

An added benefit of spin dressing is that the 3 He and neu- 
tron spins can be held closer to parallel on average, thereby 
reducing the rate of absorption of neutrons by the 3 He mag- 
netometer atoms. As discussed in Ref. [6], by modulating the 
dressing field amplitude about the point where (§„ ■ §3) = 1, 
the scintillation light will be modulated at twice the dressing 
field modulation frequency. If (S n • S3) + 1, the scintilla- 
tion light will have a frequency component at the modula- 
tion frequency. Feedback can then be employed on the DC 
component of the dressing field amplitude to drive the mod- 
ulation frequency harmonic of the scintillation light to zero 
which will hold the neutron spins parallel to the 3 He spins 
on average. This increases the effective coherence time while 
maintaining the same average signal amplitude (scintillation 
rate) compared to simple Larmor precession in a static mag- 
netic field. 

This proposed implementation of the dressed spin tech- 
nique requires that the dressing parameters be modulated to 
infer the critical dressing condition; a «EDM would be evi- 
dent by a change in the critical dressing condition as a function 
of the electric-field direction relative to the applied magnetic 
field. To date, no high-accuracy measurements or full calcula- 
tions of the influence of relatively slow modulation on the crit- 
ical dressing parameters have been performed. Furthermore, 
the theory for the effects of static field magnitude and mis- 
alignment, presented in Ref. [6] as perturbative effects, have 



not been verified. The goals of the work presented here are 
to experimentally verify the theory presented in Ref. [6] at a 
level of accuracy suitable for the «EDM experiment and to ex- 
tend the details of the calculations to encompass the possible 
effects of modulation of the dressing parameters. 

Spin dressing uses an oscillating magnetic field detuned to 
high frequency (the dressing field) to effectively change the 
gyromagnetic ratio of a spin-'/2 particle [13]. The Hamilto- 
nian for this system can be written as 

H — -yBos 7 + tiUtja ' a H — — d s x (a + a ' ) , ( 1 ) 
2 V<ra) 

where (ri) » 1 represents the average number of photons 
in the coherent, oscillating spin-dressing field, Bj is the am- 
plitude of the dressing magnetic field, Dd is the angular fre- 
quency, Bo is the static magnetic field along the z direction, 
and y is the undressed gyromagnetic ratio for the spin of in- 
terest. The operators and a correspond to the spin angular 
momentum and photon annihilation operators, respectively. 
This Hamiltonian can be rewritten by defining the quantities 
y = yBo/a>d and x = yB^/cod as in Ref. [14], making Eq. 1 
become 

H „ t „ y„ x ^ „ „ + 

a' a — — 0% h (cr-a + o~-a 

TkOd 2 4#) 

+&+a + &+a') , (2) 

where cr, are the usual Pauli spin matrices. The eigenvalues 
of this Hamiltonian, when diagonalized in the \n, ±) basis, are 
shown in Fig. 1 . The difference in energy between the two 
eigenstates around E — is then \hy'Bo\, where y' is the 
dressed gyromagnetic ratio of the spin, which depends on x 
and y. The relative sign of y' is determined by both the initial 
sign of y and by diabatically following the states as x is in- 
creased. In the limit of y <K 1, the solution yields y' = yJo(x), 
where y' and y are the dressed and undressed gyromagnetic 
ratios, respectively, and Jo is the zeroth order Bessel func- 
tion [13]. The ratio of the dressed to undressed gyromagnetic 
ratios, y' /y, is shown in Fig. 2 as a function of x for multiple 
values of y. Because of the structure of the Hamiltonian, the 
ratio ofy'/y is identical for y > and y < 0. 

In the specific case of the proposed «EDM experiment, we 
are interested in the amplitude of the dressing field that yields 
the critical dressing, at which point the two dressed gyromag- 
netic ratios are equal to each other. For y « 1, this occurs 
when y' 3 = y' n or y 3 -/ote) = y n Jo(x„), where x t = yiB d /u> d . 
Because 73 « 1.1 \y n , the point at which the dressed gyro- 
magnetic ratios are equal is when |jC3 1 ~ 1.31 (or equivalently, 
\x n \ « 1.18) in the limit wheny <*c 1. 

In this paper, we experimentally verify the spin-dressing 
predictions described above and in Ref. [6], with particular 
focus on the critical dressing point. There have been two pre- 
vious measurements of spin dressing in 3 He, one in a cell [15] 
and one in a beam experiment [14]. The latter is similar to 
this work and employed many of the same components. How- 
ever, several important changes were made that brought an 
improvement in accuracy and allow for a detailed compari- 
son between theory and experiment at a level relevant for the 
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proposed «EDM experiment described in Ref. [6]. The appa- 
ratus and the improvements are detailed in Sec. II. By numeri- 
cally solving the Hamiltonian (Eq. 2) for conditions appropri- 
ate to our experiment and comparing the results with our data, 
we verify the perturbation calculations presented in Ref. [6]. 
These results are presented in Sec. III. We also search for ef- 
fects that would prevent one from simply approximating the 
effect of modulating the amplitude of the dressing field as be- 
ing equivalent to quasi-statically changing the spin-precession 
frequency. With our current level of precision, such effects 
were not conclusively seen. Our inability to experimentally 
detect such effects in a system where the dressing parame- 
ters are varied rapidly are supported by numerical simulations 
and theoretical analysis. These calculations are discussed in 
Sec. IV. 

Because this work will focus exclusively on 3 He, we will 
henceforth use y and y' to refer to the dressed and undressed 
gyromagnetic ratios of the nuclear spin of 3 He. However, this 
spin-dressing technique can, in principle, be applied to any 
spin-'/2 particle. 



II. APPARATUS 

We investigate spin precession using Ramsey's method of 
separated oscillatory fields (SOF). A schematic drawing of 
the apparatus is shown in Fig. 3. In the main experimental 
region, a large solenoid produces a homogeneous, static mag- 
netic field in the +z direction. In addition to this solenoid, 
there are five separate magnetic field coils: two n/2 rotation 
(or "spin flip") coils, two "compensating field" coils, and a 
"dressing field" coil. The role played by each coil is described 
briefly below, and the construction and operation of each coil 
is detailed in the following several subsections. These coils 
are contained within a ferromagnetic shield to minimize the 
effect of external magnetic fields. 

A beam of spin-polarized 3 He atoms from a source (dis- 
cussed below) enters the experimental region from the left in 
Fig. 3. The first and last coils encountered by the 3 He are the 
spin flip coils, which are used to apply oscillating magnetic 
fields that induce n/2 spin rotations. The first n/2 coil tips the 
3 He spins from +z to the x-y plane. The spins then precess 
about the static magnetic field as they travel to the second 
n/2 coil where the spin component perpendicular to the ap- 
propriate rotating component of the oscillating field is flipped 
into the -z direction. The detector system, which consists of 
an analyzer magnet and residual gas analyzer (discussed be- 
low), detects only those 3 He atoms with their spins aligned 
along +z. The phase of the oscillating field in the two coils is 
the same; therefore the transmitted intensity is proportional to 
5 [l - cos(0 - </> s f)\, where (/> is the total phase accumulated 
by the spins as they precess about the +z axis in between the 
two n/2 coils, and (f> s f is the phase accumulated by the spin- 
flip field during that time. The minus sign stems from the fact 
that we detect 3 He with spin projection along +t Because 
the spin flip fields are oscillating and linearly polarized (as 
opposed to rotating or circularly polarized), our apparatus is 
completely insensitive to the sign of the gyromagnetic ratio. 



For this reason, we will take y > for clarity and note that the 
results of this experiment are independent of the sign of the 
gyromagnetic ratio. 

As can be seen from the above discussion, the minimum 
transmission occurs when the phase of the spin-flip field ((/> s f) 
equals the phase accumulated by the spins ((/>) during their 
propagation between the n/2 coils. The phase accumulated 
by the spin flip field is given by <p s f = cuL/v, where to is the 
spin-flip field frequency, L is the distance between the two 
spin-flip coils, and v is the velocity (of a particular group of 
atoms in the beam). The phase accumulated by the spins is 
given by 

rv WBW * + ?[r 2 T u--)* 

J -l/2 v 7r[J. (L/2+ i sf ) \ yj v 

Jl/2 \ 7) v 

where l s f is the effective length of the oscillating magnetic 
field generated by each of the spin flip coils and B(z) is the 
value of the static magnetic field at position z. Here, z = 
is the geometric center of the experiment. The first integral 
in Eq. 3 represents precession of the spins in between the n/2 
spin flip coils and the last two integrals in Eq. 3 represent the 
detuning of the static field from the resonance condition (a> = 
a>o = yB) while the spins are located within the spin flip coils. 

Centered in the experiment is the dressing coil, which ap- 
plies an oscillating magnetic field in the +x direction that 
serves as the dressing field. Because y' < y, the principal 
effect of the dressing field is to slow the spin precession down 
and therefore reduce cp. Because the dressing field only ex- 
tends over a small length (~ 18 cm) compared to the distance 
between the n/2 spin-flip coils (~ 51 cm), the value of y' is 
dependent on z and therefore appears as a function of z in 
Eq. 3. Midway between the dressing coil and the n/2 coils are 
two compensating coils that apply an additional, static field in 
the +z direction. This additional static field makes the spins 
precess faster while in these coils and therefore increases (p. 
One may then adjust the compensating field strength for a 
given dressing field so that the total phase difference between 
the spins and the spin flip field satisfies A(f> = (p - <f> s f = 0. 
As seen by equating <p s f = ojL/v with the expression for (p 
(Eq. 3), this condition is achieved independent of v. If A(p + 0, 
then cos 2 A(f> is not unity and interpreting the transmitted in- 
tensity as a frequency shift requires detailed knowledge of the 
beam velocity spectrum. This problem is avoided in our mea- 
surements by determining the compensating coil current that 
forces A<p = 0. 

The effects of varying the dressing field strength are deter- 
mined as follows. First, with no DC current applied to the 
compensating coils and no AC current applied to the dressing 
coil, the static field is tuned such that A(/> = for a fixed cj. 
Then, with the dressing field frequency set to a specific value 
(for a set of measurements), the dressing field amplitude is in- 
creased in discrete steps. At each step, a current is applied to 
the compensating coils to bring the 3 He transmission back to 
a minimum. Thus the n/2 coils are operated continuously at 
their proper frequency for a given static field, unlike Ref [14] 
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FIG. 3. (Color online) A schematic drawing of the apparatus used to measure the 3 He spin precession (see text for description). Not shown are 
the atomic beam source, the polarizing magnet, the analyzer magnet, and the RGA used to prepare and detect the 3 He. The coordinate system 
has the beam propagating in the z direction, the x direction pointing up, and the y direction pointing out of the page. 



where the spin-flip field's frequency is varied in an attempt to 
determine the frequency shift due to the dressing effect. 



A. Atomic Beam Source and Polarizer 

The 3 He atomic beam is produced by a 1.2 cm diameter 
tube packed with 1 mm capillaries, each 3.5 cm long. The 
nozzle is held at 1.4 K by use of a pumped liquid helium cell. 
3 He gas is delivered to the nozzle through a cooled 15 cm 
long capillary, which is in turn connected to a sequentially- 
cooled 2 mm ID tube, supplied with 3 He at room temperature 
with pressure 1.5 torr. The beam is directed through aim 
long quadrupole spin state selecting magnet (polarizer) con- 
structed from NdFeB permanent magnets (field at pole surface 
of 7.5 kG), with an effective aperture of 1.5 cm. The polarized 
atom flow rate at the output can exceed 10 14 atoms/sec. 

The polarized 3 He beam is transported through the appa- 
ratus inside a "standard length", 48" long, 32 mm diameter, 
Pyrex tube. The dressing coil and the two compensating coils 
are mounted directly onto this Pyrex tube. Inside this portion 
of the apparatus, the polarization achieved is approximately 
95%. The fact that it is less than 100% is due to the 3 He back- 
ground pressure, which is in turn a result of the inefficiency in 
pumping 3 He. 



B. Solenoid, Magnetic Shield, and the n/2 coils 

The homogeneous static magnetic field is produced by a 
solenoid that is surrounded by a 1 mm thick, Co-Netic AA 
stress annealed shield. The solenoid has 780 turns wound on 
a mean diameter of 15.36 cm and is 90.5 cm long. The shield 
has a length of 83 cm and has an inner diameter of 26.5 cm. 

The two n/2 coils were identically constructed and sup- 
ported by the G-10 glass epoxy tube onto which the solenoid 
was wound. Each is wound using 1 mm diameter wire 
wrapped on a plywood form in the shape of a 1 .9 cm by 7.6 cm 
rectangular solenoid. The coil is split axially into two identi- 
cal series driven sections separated by 32 mm, each with 40 
turns. The axis of this coil is mounted so as to produce a field 



directed perpendicular to the beam, with the short (1.9 cm) di- 
mension of the coil being aligned parallel to the beam. The 
split coil is then sandwiched between two 15 cm diameter 
disks of 0.4 mm thick Cu. A 32 mm diameter hole is cut in 
the center of each disk to allow the beam tube to pass through. 
To facilitate mounting around the beam tube, the Cu disks are 
each cut in half along a diameter and then held together with 
Cu tape. The cuts are orientated parallel to the windings; eddy 
current flow is largely parallel to the cut and the shielding is 
not affected. Mapping the oscillating field with a small pickup 
coil shows that the field maximum at the center of the coil, on 
the beam axis, falls by a factor of two at the Cu sheet, and at 
farther distances falls exponentially with characteristic length 
of order of the hole radius. Sensitive measurements with the 
pickup loop showed insignificant leakage field from the 7r/2 
coils in all regions external to the Cu shielding disks. The two 
coils are spaced by 5 1 cm. 

The coils are connected in series, which ensures equal cur- 
rent in both with no relative phase shift. In addition, a se- 
ries capacitor (determined by chosen frequency) is used to 
cancel the reactive component of the coils' impedance. The 
current is measured using a series connected 10 Q, 2 W re- 
sistor. The current is supplied by a National Semiconductor 
LM4700 30 Watt audio amplifier IC, which is driven with a 
sinusoidal voltage from an Agilent 34401 arbitrary function 
generator. The spin-flip field amplitude was adjusted to give a 
minimum transmission on resonance. Note that the transmit- 
ted 3 He beam intensity as a function of the spin-flip field am- 
plitude on resonance provides a measure of the Fourier power 
spectrum of the beam velocity distribution. This spectrum can 
be transformed to the velocity probability distribution which 
has a peak at approximately 140 m/s. The n/2 condition typi- 
cally required a current of a few hundred mA. 



C. Analyzer and Detector 

A second state selecting magnet identical to the polarizer is 
used as a spin analyzer. The input aperture of the analyzer is 
located approximately 180 cm from the output aperture of the 
polarizer. The Pyrex tube connects the polarizer and analyzer 
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vacuum systems which are otherwise separate. Guide fields 
from the polarizer (analyzer) output (input) are produced by 
two single pancake coils, each 45 cm in diameter, located mid- 
way between the output (input) apertures and the magnetic 
shield. 

The detector is a Stanford Research Systems Residual Gas 
Analyzer (RGA) model 100 with a channel electron multi- 
plier. The RGA is mounted on the side port of a conflat tee 
section, with the RGA axis perpendicular to the beam axis. 
The beam enters the tee through a 15 cm-long, 1.5 cm ID tube 
and strikes a blank flange mounted on the opposing port. In 
effect, the tee acts as a pressure buildup volume, allowing the 
3 He pressure to increase by nearly a factor of 20 relative to 
that which is obtained when the beam is delivered directly to 
the RGA. The RGA is set to mass 3, and the pressure is read 
from the RS-232 port at a maximum rate of about 8 Hz by 
use of a simple MS QUICK BASIC program. The signal-to- 
noise ratio is about 100 per reading, with a background signal 
corresponding to unpolarized 3 He or gas which reads as mass 
3. 

In order to determine the compensating coil current re- 
quired to give the minimum 3 He transmission, the drive fre- 
quency for the n/2 coils is modulated at 2 Hz with a deviation 
of 100 Hz. This is slightly smaller than the expected linewidth 
for the SOF method which is given by 



the compensating coils, i.e. 

-L/2 



„ rH1 1 v 140 m/s 

8f [Hz] = — = — = 

J 2T 2L 2 x 0.51m 



140 Hz. 



(4) 



The digital reading from the RGA is converted to a voltage 
with an 8 bit digital to analog converter and is phase-locked to 
the 2 Hz signal using a lock-in amplifier. The compensating 
coil current can then be readily set to null the detected signal 
at the modulation frequency. In the vicinity of this minimum, 
the 3 He transmission probability only depends on the phase 
difference to the second order [i.e., cos(A^) « 1 - (A0) 2 ]. 
Therefore, when the signal is set to null at the modulation fre- 
quency, it is an indication that the average phase has been reset 
to zero. For the measurements presented here, the adjustment 
of the dressing parameters and the setting of the compensating 
coil current and its measurement were done manually. 



D. Compensating Coils 

The two compensating coils are symmetrically located and 
wound directly on the beam tube and placed between the 
dressing coil and the two n/2 coils, approximately 6 cm away 
from the n/2 spin flip coils. Each coil has 1 1 turns of 1 mm 
diameter wire and an effective diameter of 1 .65 cm. Equal cur- 
rent is applied to the coils by connecting them in series, and 
their fields are orientated in the same direction. A Fluke 87 
True RMS Multimeter that was calibrated against an Agilent 
34401 A Digital Multimeter is used to measure the compensat- 
ing current, which is the current required to return the overall 
phase to the minimum transmission condition. 

The apparatus is sufficiently long that the field integral of 



XL/Z 
B c (z)dz « nnol 
L/2 



(5) 



where B c is the compensating field as a function of z, n is the 
total number of turns in the two coils, and / is the current, is 
nearly independent of the coil radius and their position in the 
apparatus. The long-distance field does penetrate into the spin 
flip and dressing coil regions and is taken into account in the 
analysis. 



E. Dressing Cosine Coil 

A cos 0-like (or cylindrical saddle-shaped) coil [16] placed 
midway between the two n/2 coils is used to apply the dress- 
ing field in the x direction. It was constructed using 1 mm 
diameter wire wound on a length of 3 in. (7.6 cm) ID PVC 
tube into which appropriate grooves had been machined on 
the outer surface. Each layer of the coil has nine axial current 
rungs per quadrant, which are then interconnected by a series 
of azimuthal current paths at either end (forming two spiral 
sets of saddle-shaped paths on the outer surface of the former). 
A total of three complete windings are employed, resulting in 
a mean coil diameter of 8.5 cm. The set of azimuthal current 
paths at either end of the coil each span an axial distance of 
1.6 cm and were positioned to give a mean length to diameter 
ratio of 2 and a total coil length of 18 cm. As designed, the 
calculated central field produced by this coil is 4.3 G/A in the 
absence of shielding [16, 17], and its measured inductance is 
0.3 mH. 

The cosine coil is driven in a fashion similar to the n/2 coils 
with the same model generator and amplifier. However, it is 
driven with much higher currents, up to 4 A RMS maximum. 
Again, a series capacitor is used to tune out the reactive com- 
ponent of the coil impedance. The current is monitored using 
8 parallel-connected 1 1 Q. 2 watt carbon composition resis- 
tors. The resistance was calibrated as a function of current 
and frequency to check for possible nonlinear heating or other 
effects. There is no frequency effect, however a small heating 
effect with increasing current was taken into account in the 
analysis. An Agilent 34401A True RMS Multimeter is used 
to measure the voltage across the current monitor resistor for 
the dressing measurements and calibrations. 

Provisions were made to add a DC current in addition to the 
oscillating current in the dressing coil, allowing both the de- 
termination of the absolute average angle between the dress- 
ing field and static field, and to study the effects of intentional 
misalignment. This was accomplished by use of a choke (in- 
ductor) of approximately 30 mH, connected to the node be- 
tween the capacitor and dressing coil. The choke presents a 
reactive impedance of 500 Q at the lowest dressing frequen- 
cies, much higher than the dressing coil impedance. Only the 
current through the dressing coil was sent through the monitor 
resistor. There was no net effect on the AC current due to the 
presence of the choke. 
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III. SPIN DRESSING RESULTS 

Data were acquired at undressed Larmor precession fre- 
quencies of a>o/2n = 3.99 kHz and o>q/2jt = 9.36 kHz. Fig- 
ure 4 shows the compensating field required to maintain min- 
imum transmission as a function of dressing parameter x. 
The compensating field strength is expressed as its average 
strength over the entire length of the apparatus divided by the 
magnetic field required for the undressed Larmor precession 
frequency, i.e. 

Bcav 1 C L/2 

-F = 1T7 ^c(z) dz , (6) 

-do dqL J-l/2 

where Bo = <*>qIj- For each undressed Larmor precession fre- 
quency, multiple sweeps of the intensity of the dressing field 
with different frequencies were taken, allowing the phase shift 
to be measured as a function of both dimensionless parame- 
ters x and y. Note that the functional form of the prediction 
mimics the behavior shown in Fig. 2, but is inverted. 

In the region between the coils, we must know the field pro- 
files accurately in order to make a prediction. Computation of 
the field profiles due to the solenoid and compensating coils 
was made using a finite-element solver in the presence of the 
magnetic shielding. These calculations were verified using 
data where, without dressing, the magnitude of the solenoid 
field was reduced to compensate for the additional field ap- 
plied by the compensating coils to maintain minimum trans- 
mission, paying special attention to the last two integrals in 
Eq. 3. The dressing coil field profile was mapped with a pair of 
orthogonal 7 mm diameter (axial and transverse) search coils 
and with a Hall probe, both with and without the magnetic 
shielding, and as a function of frequency. Not only does the 
presence of the magnetic shield change the field calibration, 
but the calibration is dependent upon the frequency of the field 
generated by the coil. At these relatively low frequencies, the 
skin depth is expected to be much larger than the thickness of 
the surrounding magnetic shield. However, because of vari- 
ous geometric effects, the penetration is not merely a function 
of the skin depth alone, resulting in a more complicated de- 
pendence on frequency [18, 19]. The field profiles used in this 
analysis are shown in Fig. 5, along with the measured depen- 
dence of the dressing coil field on frequency. 

In addition to incorporating the influence of the shield on 
the dressing coil calibration, we include the possibility that the 
dressing field is not perpendicular to the field of the solenoid 
and compensating coils. By applying a static field with the 
dressing coil rather than an oscillating field, the fields due to 
the solenoid and the dressing coil add in quadrature. Let us 
define the x direction as being parallel to the dressing field. 
The first integral in Eq. 3 then has y'(z) = y (because there is 
no dressing in this case) and 

B(z) = V(#o sin 9 + B c ) 2 + (B ds + B Q cos 0) 2 , (7) 

where B ds is the static field applied by the dressing coil, B Q 
is the field applied by the solenoid, B c is the field due to the 
compensating coils, and 6 is the angle between the solenoid 



field and the dressing coil field. By recording the compensat- 
ing coil current required to keep the transmission minimum at 
multiple DC dressing coil currents, we can determine the best 
fit angle between the dressing coil field and the solenoid field, 
which is 88.6 +0.1 deg. As shown in Ref. [6], the primary 
effect of this misalignment when the dressing is active is to 
only dress the component of the spin precession along the axis 
perpendicular to the dressing field, while the spin continues to 
precess normally about the component of the static field paral- 
lel to the dressing field. This effect is included in this analysis. 
However, further effects not considered in Ref. [6] can result 
from misalignment and these are discussed in Sec. IV. 

For the data shown in Fig. 4, we diagonalize the Hamil- 
tonian of Eq. 2 with the misalignment discussed above and 
determine the effective gyromagnetic ratio by taking the dif- 
ference of the two central energy eigenstates and dividing by 
HBq. This effective gyromagnetic ratio is then determined at 
each position z in the apparatus. By assuming that the spin- 
precession frequency changes quasi-statically, this dressed gy- 
romagnetic ratio appears as a function of z, i.e. y'(z) in Eq. 3. 
In order to simplify the integrals, we make the assumption that 
in the region where the correction field B c (z) is non-negligible 
[i.e., B c (z)/Bq » 1], the precession frequency is set by the 
undressed y, i.e., y'(z) ~ y. We note that this occurs when 
B c ,av ~ 0.1 B , because B c (z) has a much smaller spatial ex- 
tent then Bq. In these regions, the fringe field produced by the 
dressing coil is < 1% of the maximum (dressing) field and so 
the maximum local value of x incurred is only x ~ 0.08 when 
lo q /2tt = 3.99 kHz and u d /2n = 5.12 kHz. 

Using these assumptions, we calculate the required com- 
pensating field to maintain minimum transmission. These the- 
oretical curves are shown in Fig. 4. Note that there are no 
tunable parameters in the curves shown; all unknown param- 
eters in the theory are determined by calibration. Discrep- 
ancies between the data and the theory are less than 4% for 
1 < x < 6 when the undressed Larmor precession frequency 
is a>o/27T = 3.99 kHz and less than 2% for all values of x mea- 
sured when cl>q/27t = 9.36 kHz. The residuals are not Gaus- 
sian and do show some systematic error, which is most likely 
caused by inaccurate mapping of the magnetic fields. Any 
greater precision would require maps of the magnetic fields to 
be more accurate than 1%. 



IV. BREAKDOWN OF THE QUASI-STATIC 
APPROXIMATION 

The primary assumption made in the analysis presented 
above is that the spin-precession frequency changes quasi- 
statically. One particular region where this assumption can 
break down is where y'(z) — > 0. In such regions, if the align- 
ment between the dressing field and the static magnetic field 
is not perfect, an avoided crossing of the eigenvalues of the 
spin dressing Hamiltonian appears, as illustrated in the inset 
of Fig. 1 . The size of this avoided crossing corresponds to un- 
dressed precession of the spin about the static magnetic field 
that is parallel to the dressing field. 

In the sudden, semi-classical approximation used above, the 
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FIG. 4. (Color online) Left: Average required (integrated) compensating field strength in units of the undressed magnetic field for an undressed 
precession frequency of 3.99 kHz and dressing field frequencies 5.12 kHz (y = 0.779, squares), 8.59 kHz (y = 0.464, circles), 16.4 kHz 
(y = 0.243, upward triangles), and 35.8 kHz (y = 0.1 1 1, downward triangles). Right: Average required compensating field strength in units of 
the undressed magnetic field for an undressed precession frequency of 9.36 kHz and dressing field frequencies 13.75 kHz (y = 0.681, squares), 
16.44 kHz (y = 0.569, circles), and 20.50 kHz (y = 0.457, upward triangles) 
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FIG. 5. (Color online) (Top) Calculated magnetic field profiles for 
the solenoid coil (short-dashed, green curve) and compensating coil 
(long-dashed, blue curve), along with the experimentally determined 
profile for the dressing coil in the presence of the magnetic shield- 
ing (solid, red curve). (Bottom) The magnitude of the dressing field 
vs. frequency of the oscillating current in the presence of magnetic 
shielding. 



spin precesses about the vector sum of the static field point- 
ing in the same direction as the dressing field (B x , in the x 
direction) with angular frequency yB x and the component of 
the static magnetic field orthogonal to the axis of the dressing 
field {B z , in the z direction) with angular frequency y'(z)B z . 
Thus, in this approximation, the phase accumulated through 
the region of interest is 



A0 




B z vt\ +(yB x ) 2 dt, 



(8) 



where y is the undressed gyromagnetic ratio, y'(z) is the 
dressed gyromagnetic ratio, v is the velocity, and t — is the 
time at which y'{z) = 0. Here, the first term represents the pre- 
cession about z, where the instantaneous frequency is given by 
w z (?) = j'(z(t))B z = c ^-^-B z vt, and the second term represents 
the instantaneous precession frequency about i. When treat- 
ing the full precession frequency as merely the (frequency) 
difference between the energy eigenvalues of the of two states 
as shown in Fig. 1, we implicitly make the assumption that 
(jj(-6f) - -coiSf), where 6t is an infinitesimal time step from 
t — and u is the full angular precession frequency. Semi- 
classically, this amounts to time-reversing the spin dynamics 
regardless of the magnitude of the orthogonal component B x . 
Quantum-mechanically, we assume that we are in the strong 
diabatic transfer limit, because in this picture a state |J,) re- 
mains in |X) regardless of the energy of that state as a function 
of time. 

Rigorous treatment of the effect of B x complicates the prob- 
lem significantly for two reasons. First, because B x is along 
x, the spin will precess outside of the x-y plane in the region 
where y'B z <K yB x . After traversing this region of interest, 
the spin can be left in a state with a non-zero value for (S z ). 
Second, if the spin is moving sufficiently slowly through the 
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FIG. 6. (Color online) Results of the simulation of the approximate 
Hamiltonian (Eq. 9, see text). In the top pane is the final expecta- 
tion value for S z (— V2 S (S z ) < V2) an d the bottom pane contains 
the total phase swept out as the spin traverses the avoided crossing. 
Note that the simulation is highly dependent on the initial state. The 
solid, red lines correspond to an initial state [Eq. 10] with = 0, the 
blue, dashed lines correspond to <f> = n/2, and the dotted, green lines 
correspond to <p = n/4. 



region of interest, the spin-precession vector will adiabatically 
follow the vector sum of yB x and y'B z . To be specific, at the 
start of the region of interest, when y'B z » yB x , the spin- 
precession vector will point along +2. As y' — > 0, the spin- 
precession vector will point along x, as it will only be pre- 
cessing about yB x . As y' becomes negative and its magnitude 
grows, the spin-precession will realign with the -z axis. The 
spin therefore accumulates an additional phase of ji radians as 
it started along +z and ended along -f . This additional phase 
can be approximated by having the spin rotate around a pseu- 
domagnetic field pointing in the y direction and appears to be 
similar to a geometric phase effect. 

Quantum mechanically, this problem is similar to a Landau- 
Zener tunneling process, in which the Hamiltonian for the re- 
gion of interest can be well approximated by 



H = 



1 



'-%B z vt yB x \ 
yB x ^B z vt) 



(9) 



In our case, the initial state lies somewhere in the x-y plane, 
e.g. 



-10/2 



IT> + ^ /2 ||>) , 



where <p represents the azimuthal angle. The state is then 
evolved through the avoided crossing and the additional phase 
accrued is examined. To the authors' knowledge, no analytic 
solution exists to the problem with this initial state. Nonethe- 
less, we expect that the figure of merit remains the same as for 
the standard Landau-Zener problem, namely 



y 2 B 2 « — B z v . 
dz 



(11) 



in order for the spin to proceed in the diabatic limit and accu- 
mulate no additional phase through the region where y' — > 0. 
In the present case, B z = Bo sin 9 and B x - Bo cos 9, where 
9 = 88.6°. Given the angle, we further approximate B z * B . 
The static field Bo can be determined through the undressed 
precession frequency, e.g. Bo = u>o/y = fo/(y/2n) « 
10kHz/(3.243 kHz/G) * 3 G. As discussed in Sec. IIB, the 
average velocity for the 3 He was measured to be v = 140 m/s. 
For the purposes of this estimate, let us consider the case 
where the amplitude of the dressing field is large enough such 
that the dressing parameter at z — (the maximum value of 
x(z)) is given by x(z = 0) > 2.5. Using the amplitude profile 
shown in Fig. 5, the 3 He dressed gyromagnetic ratio will be- 
come zero at approximately z ~ +5 cm. The rate of change of 
y' around these points can also be estimated from the ampli- 
tude profile in Fig. 5. With d(y'/2n)/dz ~ 25 (kHz/G)/m, 

Idy'Y 1 y 2 B Q cos 2 9 , 

-7- — ~3.5xl(T 2 ; (12) 

\dz j v 

so we are indeed in the diabatic limit. 

Simulations were made of the full time-dependent 
Schrodinger equation with the Hamiltonian of Eq. 9 to 
quantify the effect of this misalignment. The initial state 
(Eq. 10) is propagated from the dimensionless time "f = 
(dy' /dz)B z v) l/2 =-100 to?' = 100, and the total accumulated 
phase and final values of (S z ) are computed as a function of 
yB x . These results are shown in Fig. 6 for various initial states. 
The final state depends strongly on the initial state. However, 
for the diabatic case, we empirically find from simulation that 
an order of magnitude approximation of the effect is given by 



(dy'Y 1 y 2 B cos 2 6 



(13) 



(10) 



which implies an additional shift of approximately 0.11 rad 
as the spin traverses the region of interest when x(z = 0) > 
2.5. This phase shift is less than 0.1% of the total phase (f> 
accumulated by the spin as it traverses between the spin flip 
coils with no dressing applied. 

Even when the dressing field is perfectly orthogonal to the 
static field, an additional phase shift to the quasi-static ap- 
proximation can arise. As the spin enters the dressing coil, 
it is subject to a changing magnetic field due to the changing 
intensity profile of the dressing field (see Fig. 5). This effect 
is well known in interferometer experiments, and a similar ef- 
fect was discovered by Millman [20]. In the rest frame of the 
atom, the 3 He experiences an amplitude-modulated dressing 
field that changes the dressing field's effective spectral pro- 
file causing dressing at multiple frequencies, an effect that is 
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not included in the quasi-static approximation. In order for 
the influence of this modulation to be small, the amplitude of 
the dressing field must change more slowly than the overall 
precession, so that the spin accurately follows the field as re- 
quired by the quasi-static approximation. Simulation of the 
full Schrodinger equation for our field profiles indicates that 
for x < 3, this Millman-type effect scales roughly as 



\or d ) az 

Note that this effect will become stronger as the velocity is 
increased, in contrast to the situation for the avoided cross- 
ing effect discussed above. For this reason, we expect that 
this Millman-type effect will be the predominant deviation 
from the quasi-static behavior encountered in our experiment. 
One of the largest phase shifts predicted for this experiment 
is A(f> * -0.8 rad at w /2tt = 3.99 kHz, w d /2;r = 5.12 kHz, 
x(z — 0) = 3, and v = 140 m/s. This phase shift is slightly less 
than 1% of the total phase cp, on the threshold of what might 
be detectable in our data. 

Searching for these effects in our experimental data requires 
a significant computational effort. Because these are purely 
dynamical effects and no analytical solution exists, numerical 
integration of the Schrodinger equation must be made which 
presents several complications. First, the field profiles shown 
in Fig. 5 must be used. Second, because there is a distribution 
of velocities and each of the above effects depends strongly on 
velocity, each velocity class must be simulated separately and 
subsequently averaged. Third, when y ~ 1, the phase of the 
spin dressing field will affect the final state; therefore averag- 
ing over the phase of the dressing field is also required. We ran 
simulations of our Ramsey experiment which calculated the 
transmitted 3 He intensity as a function of velocity, the dress- 
ing field phase, and x for a given value of y. We first averaged 
over the phase of the dressing field, assuming that all possible 
phases are equally likely for each velocity class. We then av- 
eraged over the velocity classes using a Maxwell-Boltzmann 
distribution with (v) = 140 m/s, which corresponds to a tem- 
perature of approximately 2.9 K. For each value of x, the sim- 
ulation tuned the strength of the compensating field to mini- 
mize the total transmission. 

Typical results of a simulation are shown in Fig. 7. In gen- 
eral, the full quantum simulation better replicates the data 
for x < 2.3, but deviates from the data for larger x val- 
ues. This deviation appears to be due to the introduction of 
the avoided crossing effect, which causes rapid oscillations 
in the transmission vs. velocity curve. These oscillations 
could easily lead to numerical errors on the 1% level. The 
additional phase shift that is accumulated in these simulations 
compared to the quasi-static approximation is due primarily to 
the Millman-type effect, as the additional phase shift roughly 
follows Eq. 14. A comparison of the sum of the squared resid- 
uals for the two predictions is inconclusive, as no theory pro- 
duces a better curve for all data sets. However, there may be 
indications that the Millman effect is still seen in the data. For 
x < 2.3 and y ~ 1, the residuals for the full simulation show 
little or no x dependence whereas the quasi-static approxima- 
tion tends to show a linear dependence of the residuals on x. 
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FIG. 7. (Color online) Results of a full simulation integrating 
the Schrodinger equation for a>(,/2n = 9.36 kHz and u> ( il2n = 
13.75 kHz. Top panel: Data, shown as points, and comparison the- 
ory. The solid blue line shows the quasi-static approximation while 
the dashed green line shows the full simulation of the experiment. 
Bottom panel: Difference between data and theory for the quasi- 
static approximation (blue, upward triangles) and the full simulation 
(green, downward triangles). 



Moreover, the sum of the squared residuals for x < 2.3 tends 
to be a factor of 3-4 smaller than that of the quasi-static case. 

While easily distinguishable, the absolute difference be- 
tween the two theoretical results is typically less than 2%. 
In order to conclusively demonstrate these effects experimen- 
tally, calibration errors must be suppressed to within 2%. Sev- 
eral possible systematics, such as external magnetic fields 
leaking through the magnetic shielding, improper determina- 
tion of the magnetic field profiles, and other such effects can 
easily account for these discrepancies. However, absolute de- 
termination of these systematics is beyond the scope of this 
work, and therefore conclusions regarding these deviations 
from the quasi-static condition cannot be made. 



V. CONCLUSIONS 

In this work, we verified the spin-dressing Hamiltonian 
(Eq. 2) to a level of 2% in the region of interest for a pro- 
posed neutron EDM experiment, verifying the perturbative 
corrections contained in Ref. [6] to the same level. We also 
performed simulations to search for breakdown of the quasi- 
static approximation used above and in [6]. Ultimately, fur- 
ther experimental work (in which systematics are controlled 
to an absolute level of 1% or better) would be required to 
definitively identify effects associated with the breakdown of 
this approximation. Alternatively, it could be studied by con- 
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trolling the atomic velocity distribution. In particular, faster 
atoms are expected to enhance the magnitude of departures 
from quasi-static behavior. This approach is not feasible us- 
ing our present apparatus; in effect we are unable to polarize 
atoms traveling much faster than those used in this work be- 
cause of limitations associated with the magnetic fields in the 
polarizer and analyzer assemblies. One might be able to study 
the breakdown of the quasi-static approximation using a cell- 
type experiment, such as that presented in Ref. [15]. In this 
case one would use a rapidly modulated dressing field to ex- 
plore appropriate regions of parameter space. 

An understanding of the implications of the breakdown of 
the quasi-static approximation is required for the proposed 
neutron EDM experiment [6]. In order to study the effect on 
the proposed experiment, we simulated the effect of modu- 
lation of the spin dressing using realistic parameters. In the 
simulation, an undressed Larmor precession rate of 20 Hz 
was dressed with a modulated dressing parameter, x = x c + 
x m cos(w m f), where x c = 1.31, x m = 0.1, and (o m /2n = 1 Hz. 
The dressing frequency was set to 2 kHz. Note that this places 
us in the regime where y« 1, which greatly suppresses the 
Millman-type effect discussed in Sec. IV. Also, because the 



region near x = 2.3 is avoided, any misalignment effects are 
also suppressed. Because of these suppressions, the quasi- 
static approximation works quite well to within the numerical 
errors of the simulation. No phase lags or unexpected har- 
monics of the modulation frequency were seen in the simula- 
tion. However, more simulations or experimental work must 
be done for the final parameters used in the experiment. This 
would include simulating a square wave modulation, which 
would provide a larger signal to noise ratio for the proposed 
neutron EDM experiment. 
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